Gamma rays and neutrinos from a powerful cosmic accelerator 



Abstract 



K. Mannheim 1,2 
1 Max- Planck- Institut fur Radioastronomie 
1 Auf dem Hiigel 69, W-5300 Bonn 1 

0^ ' E-mail: p036kma@mpifr-bonn.mpg.de (INTERNET) 

G ■ 2 Universitats-Sternwarte 

Geismarlandstr. 11, D-37083 Gottingen 
E-mail: KMANNHE@gwdgvl.dnet.gwdg.de (INTERNET) 

■ Federal Rebublic of Germany 

in ; 

(February 5, 2008) 

O 

CO 

Pur 

2 ' Possibly, the powerful radio quasar 3C273 will reveal its nature as an efficient 

proton accelerator up to energies of order 10 11 GeV in the near future. It is 

■ shown in this paper that the shock accelerated protons expected to be present 
in the quasar's plasma jet induce an unsaturated synchrotron cascade with 
electromagnetic radiation emerging in the X-ray and gamma ray range. While 
(including the synchrotron emission from the accelerated primary electrons) 
the broadband nonthermal emission from 3C273 can be explained over the 
observed 18 orders of magnitude, a flattening of the spectrum at the highest 
observed energies (a few GeV) is predicted that could be falsified by the En- 
ergetic Gamma Ray Experiment Telescope (EGRET) on board the Compton 
Gamma Ray Observatory (GRO). Above ~ 100 GeV the cascade spectrum 
dramatically steepens again due to the absorption of the gamma ray photons 
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by the host galaxy's strong infrared photon field from extended dust clouds 

- in accordance with the non-detection of 3C273 by Cherenkov telescopes. 

However, neutrinos from the hadronic interactions initiating the cascade are 

not damped and reach terrestric experiments without any modification of 

their injected flux. In contrast to the neutrino flux from pp-inter actions, 

which are energetically unimportant in jets, p7-interactions generate a flat 

neutrino flux. Therefore it is emphasised that one must not simply normalize 

the expected neutrino flux by the observed gamma ray flux. Hence it is shown 

that the expected neutrino flux in the energy range relevant for underwater 

or underice detectors is much lower than assumed by many authors. On the 

other hand, with an increasing number of cosmic gamma ray sources at known 

positions, their neutrino detection should be feasible when it is realized that 

angular resolution is the crucial design property for neutrino detectors. 
95.30, 95.70.R, 98.70.S, 96.40.P 
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I. INTRODUCTION 



Recent gamma ray observations by GRO may be the clue to the long-standing problem of 
the origin of UHE-cosmic ray protons. GRO has discovered that many flat-spectrum radio 
quasars are strong gamma ray sources with energy spectral index s = log [dN 1 /dE}/d\og [E] 
above 100 MeV ranging between roughly —1.5 and —2.5. As pointed out below the emission 
indicates the presence of very energetic baryons that interact with the dense photon atmo- 
sphere in the plasma outflow emerging from the active galactic nucleus (AGN). A prominent 
example is the nearby quasar 3C 273 with index — 2.39±0.13 |T| which is investigated in this 
paper because the steep gamma ray spectrum observed may help to discriminate between 
the proton powered and other - purely leptonic - models. 

A. The nonthermal "mini-blazar" 

3C 273 is a well-known superluminal radio quasar at redshift z = 0.158 (e.g. Ref. 
0). VLBI experiments showed that radio knots propagate along a strongly collimated jet 
emerging from a stationary core with apparent opening angle $ Q b = 4°. The motion of 
the knots with respect to the core is superluminal with apparent speed /3 D b = 7/i _1 over a 
projected distance of at least r Q b = 45/zr 1 pc at 5 GHz where h = H o /100 Mpc km -1 s -1 . 
The total flux of the jet comprises the flux of several knots with increasing self-absorption 
frequency towards the base of the jet and is flat (a ~ 0) up to a spectral break in the mm 
range above which the spectrum steepens (a ~ 1). Detailed studies [§] of the variability and 
polarization of 3C 273 have disentangled the nonthermal emission component associated to 
the compact jet (the low-entropy "miniblazar" ) from the thermally emitting (high-entropy) 
sites, such as the UV and soft X-ray excess emitting "cold" gas or the infrared emitting 
dust. 

The miniblazar can be attributed to a relativistic plasma outflow from the AGN with 
embedded magnetic fields and shocks |4j . The shocks accelerate particles and the accelerated 
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electrons loose their energy rapidly due to observable radio-optical synchrotron emission 
in the magnetic field. Interferometric radio maps rather suggest a helix pattern for the 
trajectories of the radio knots and, indeed, the situation can be pictured as a topological 
equivalent to the curved shocks in the solar wind following the Parker spiral of the magnetic 
field - with the equatorial plane of the sun being bent onto a cone, so that the spiral becomes 
a helix, and the shock waves being excited by mass ejections from the accreting gas at the 
base of the jet. Remarkably, the accelerated protons and nuclei || were long thought to be 
negligible as radiative agents, since they have much lower energy losses at the same energy 
as electrons. However, acceleration thus operates longer for the heavy particles until their 
energy losses also become large. The energy lost by the energetic baryons then emerges as 
cascade radiation in the X-gamma ray band. 

Since the shock speed does not reflect the fluid speed (in which the scattering centers 
move with the Alfven speed), the speeds of knots may be different, i.e. a strong recon- 
finement shock || after an initial passage of a thick cloud may be a standing or a weakly- 
relativistic shock wave. As a consequence the emission from this shock, particularly the 
infrared/X-ray emission, could be less beamed than the radio emission, since radio emission 
is self-absorbed within the inner parsec of a jet and thus comes from further out. Beaming of 
the emitted radiation pattern results from the anisotropy arising from the Lorentz transfor- 
mation between the comoving fluid frame (with isotropic particle and photon distributions) 
and the observer's frame ("Doppler-boosting"). However, for the sake of simplicity it is 
assumed in this paper that the Doppler factor is the same along the jet. 

B. UHE protons 

The crucial question is, of course, wether protons can really reach the extremely high 
energies needed to generate a luminosity of secondaries competitive with inverse-Compton or 
synchrotron-self-Compton emission from the accelerated electrons and to be able to provide 
enough energy in the center-of-momentum frame of the colliding protons and photons to 
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create secondary particles on the mass shell. 

The maximum possible energy is readily obtained from 

^p, cool (-^p, max) ^acc (-Sp,max) (I'-O 

provided that the limit from drift across the shock with radius r±_ 

£ P ,max « e(v x B)r/c ~ 7.8 • 10 2 °A (—) eV (1.2) 



V G/ \pc / 

is not exceeded. Here t\ denotes the time scale of process i, B the magnetic field 
and P s the speed of the shock. Due to energy losses the proton distribution steepens from 
dN/dE oc E~ 2 (or even flatter, cf. Ref. [|]) to dN/dE oc E~ 3 when 

^PjCooK-^p,!}) ^e.cooK^Cjb) ^exp (^•^) 

where t exp denotes the dynamical time scale of the expanding jet. When E Vt \> > E PtTnaiX , 
no steepening occurs and the proton distribution turns over with a Bessel-function type 
behaviour at -E P)ma x- Since 

L p _ u p t e>coo i ^ Min [E Pih , -E p , max ] , 

T n 4- TP ' 

L/ e U e tp )C ool Pjb 

a small E P)iaax costs a great value of rj = u p /u c to obtain comparable luminosities from 
protons and electrons, respectively. Note that rj in Eq. (|1.4j) refers to the energy density 
ratio at extremely relativistic energies, so that the Galactic value rj(> GeV) ~ 100 indeed 
lets one expect the possibility that L 7 ^> L ir . With B oc r _1 along the jet of length r 
and r small enough, so that E p ^ < E Pjmax , it follows that L p /L e oc rj, whereas further 
out L p /L e -C r). Therefore, one expects to observe gamma rays mainly from the compact 
regions of radio jets close to the core. Additionally, the emission from the nuclear jet (the 
blazar) is concentrated in a narrow lighthouse beam due to relativistic bulk motion of the 
plasma amplifying the flux at small angles to the line of sight. Thus it can be explained, 
why flat-spectrum radio sources (indicating jet sources at small angles to the l.o.s.) are the 
gamma ray sources detected in flux-limited experiments. 



Comparing 3C 279 and 3C 273 one can see that in spite of a sim ilar inferred radiation 
compactness, the maximum proton energy in the jet of 3C 273 is smaller because of the 
presence of the strong external radiation field seen as the big blue bump. For the maximum 
dimensionless energy it follows from Eq. (|1.3|) that 



1 + (<Xp 7 /er p ,syn)as 

with the target photon/magnetic energy density ratio a s = u^/ub and the factor 
(<j p7 /(Xp ]Syn ) ~ 240 accounting for the relative importance of photoproduction and syn- 
chrotron energy losses, respectively. For the intrinsic photons a Si i n ~ /5j7j$/(l < 1, 
whereas the external photons with luminosity L uv ~ 10 47 erg/s enhance the local pho- 
ton density such that a Sjex = 2L uv /(r^ + r c ) 2 B^c ~ 0.1 assuming = r c for the length 
rb of the radiative jet and the distance r c of the base of the jet to the location of the 
thermal photons (i.e. collimation out to the parsec scale). The value thus obtained 
for the proton break energy 7 P b = 3 ■ 10 10 (Tab. |) is in modest agreement with the 
value for 7 P)max obtained from Eq. fll.lp . Thus, head-on collisions with a (infrared) pho- 
ton of energy e t = 145 MeV/27 Pj b — 2 • 10 -3 eV can excite the A resonance. Note 
that the observed energies are cosmologically redshifted and Doppler blueshifted yielding 
-Ep,max(cbs) = _Dj7p max /(l + z) ~ 2 • 10 11 GeV, if there were no energy loss during propa- 
gation through the cosmic microwave background. 

Since the pair absorption cross section peaks for head-on collisions, but the external 
photons appear highly anisotropic in the comoving fluid frame, the pair creation opacity 
(oc compactness) in the very forward direction remains unaffected by the additional target 
photons. In contrast, the Thomson scattered luminosity is found to be Lt ~ a s ex L s , where L s 
denotes the apparent soft synchrotron luminosity of the jet, since the Thomson scattering is 
isotropic in the comoving frame. In models where the particle acceleration takes place much 
closer to the AGN, i.e. r c C 1 pc, the thermal seed photon energy density could exceed the 
magnetic field energy density in the jet, so that, indeed, the Compton-upscattered thermal 
photon flux can dominate over the emission from the jet itself (cf. Ref. fTTJ or for another 
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version of the same model see Ref. pj|). Very close to the AGN the extremely strong 
Compton cooling is faster than diffusive acceleration processes which are rather slow for 
reasonable diffusion coefficients. Only explosive acceleration seems possible in this case. 
On the other hand the well-ordered structure of the jet on much larger scales must not be 
destroyed by such processes, nor must there be too efficient dissipation, for otherwise the 
luminosity of the large scale jet could not be as large as observed. 

In the following Sections it is shown how the cascade respondes to a varying proton max- 
imum energy and how the infrared photons from dust clouds modify the emerging spectrum. 
Finally, the neutrino flux from 3C273 is calculated. 



II. THE UNSATURATED SYNCHROTRON CASCADE AND THE FEW MEV 

BUMP 

It is instructive to calculate step by step the development of the unsaturated synchrotron 
cascade induced by the UHE protons. For each step conservation of the total emitted power 
L oc x 2 dN 1 /dx determines the normalization of the stationa ry photon spectrum dN^/dx. 
The result of a numerical calculation is shown in Fig. [I] for the physical conditions listed in 
Tab. |I| and the detailed integral equation solved is discussed in [|l^ . 

We start with protons of Lorentz factor 7 p b = 3 ■ 10 10 producing pions of Lorentz factor 
j n = (nip/m^) 7 Pl b/5. The neutral pions decay yielding two gamma rays of dimensionless 
energy 

x = ^~ — ^7 P ,b - 180 7p , b (2.1) 
m e c 2 m e 10 

in the comoving frame (multiplication with the redshifted Doppler factor -Dj/(1 + z) yields 
the observed frequencies or energies). The local emissivity of gamma rays is Q 7 oc x _1 
and the stationary photon spectrum steepens to dN^/dx oc x~ 2 because of the 77 pair 
production losses. Here we take dN p /dj p oc 7~ 2 (shock acceleration, neglibible losses below 
7 P) b) and dN lttaxget /dx oc x~ 2 (as observed between 10 11 Hz and 10 15 Hz corresponding 
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to proton energies at threshold between 71 = m 7r /(2m c xi) = 10 12 and 72 = 10 8 , resp.). 
The reason for the flat injected flux of secondary gamma rays is the following: Compared 
to cooling on monoenergetic target photons with energy x m , which yield Q 1 oc x~ 2 above 
7 m = m 7r /(2m e x m ), cooling in the inverse power law photon distribution with iV 7itai . g0 t(> 
x) oc x" 1 yields additional target photons (x < x m ) for 7 p > 7 m and hence a flatter emissivity 
<5 7 oc x~ 2+1 |13|]. Now, we forget about charged pions, direct pair production and proton 



synchrotron emission, because the pairs from charged pion decay simply add more flux to 
the neutral pion decay cascade, direct pairs radiate mostly around TeV where - as will 
be shown- the emission is absorbed by dust infrared emission and the proton synchrotron 
emission contributes only for a s < 0.004. 

The first synchrotron photon generation is now radiated by the pairs produced from the 
original neutral pion decay quanta via 

7 + 7t -> e+ + e ~ 

where 7 t denotes a soft photon at threshold energy xx t = 1 (head-on collision). The 
maximum electron Lorentz factor is given by 71 = x Q /2 and the characteristic synchrotron 
frequency is then 

Xl = 3 • 10~ U B^ 2 = 2.2 • lO 11 ^ (2.2) 

The stationary electron distribution being subjected to synchroton losses is dN e d^ e oc 7 -2 
(for 7 < 71 ) and therefore the synchrotron emissivity is Q y oc x~ 3 / 2 (for x < X\). As long as 
we are in the optical thick range, the stationary photon spectrum is steeper by one power, 
i.e. dN-y/dx oc x~ 5 / 2 . The pairs thus produced have Lorentz factors 72 = xi/2 and radiate 
the second synchrotron generation 

X2 = 3 ■ 10" 14 5 ± 7 2 = 3.6 ■ 10 8 fii (2.3) 

with stationary photon spectrum dN^/dx oc x _11//4 . However, part of this generation is 
already at energies below 
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where the jet becomes optically thin (compactness in the comoving frame I — 2 • 1CT 6 for 
3C 273. cf. Tab.l), and therefore has a spectrum dN^/dx oc x~ 7 ^. The energy x** where 
the pairs produced at x* radiate is given by 



x** = 3 • lO -1 ^ 



2 

14; 



X 

T 



0.2B ± (2.5) 



The third generation has the characteristic energy 

x 3 = 3 • 10' 14 B ±1 l = 10 3 B 7 ± (2.6) 

Now, here comes the most important point: There are two special energies involved. The 
energy x* defined by r 77 (x*) = 1 and the energy x 2 , which may be either <, = or > than 
x*. In the first case the second cascade generation is optically thin with respect to pair 
creation (no further reprocessing) and has a spectrum with index —7/4. In the second and 
in the third case, there will be a third cascade generation. The part of the spectrum from 
the optically thick range x 2 > x > x* (which is monoenergetic for x 2 — x*) has the form 
dN^/dx oc x" 15/8 over an energy band A log [re] = 2 log [x 2 /x*\ (a la 3C 279) from x** to 
x 3 and the part from the optically thin range has the spectrum dN^/dx oc x~ u / 8 emerging 
below x**. Thus, a 7 ~ 0.9 and ax — 0.4 ignoring that the superposition of individual 
cascade generations somewhat smears out these values. 

The case x* x 2) relevant for 3C 273 where x 2 follows the reduction of 7 Pi t>, leads to a gap 
between the maximum of the third generation at x** x% and the maximum of the second 
(partially absorbed) generation at x*. It is not possible to predict the accurate shape of this 
gap, because it depends sensitively on the upper turnover of the injection modulo all decay 
kinematics and successive synchrotron smearing (x oc 7 2 ). The numerical result (assuming 
B± ~ B h so that DjX 3 (l + z) ~ 2, that is 1 MeV) with injection of an exponential turnover 
is shown in Fig.l for 3C 273 yielding an effective spectral index close to the observed —2.39. 
However, at a few GeV the spectrum of the decadic power rises again due to the onset of 
the second cascade generation with photon index —7/4 up to 10 TeV^ D^x* /(l + z). 



The energy X3 oc B]_ is extremely dependent upon the magnetic field strength in contrast 
to x** oc B±. Therefore one may expect a wide range of spectra in the MeV-GeV range. 
There is no way, however, to understand a few MeV bump without an accompanied intrinsic 
bump at 1-10 TeV, since the cascade generation preceding the one that makes the MeV 
bump must show up at D$x*/(1 + z). 

III. INFRARED ABSORPTION OF THE FEW TEV BUMP 

This brings us to the final argument, viz. the absorption of gamma rays outside of the 
jet by photons from the extended regions of 3C 273. 

The optical depth for scattering on a monoenergetic photon target is approximately given 

by 

r 77 (x) = r ext N 1 (x )a(x) (3.1) 

where x = hv /m e c 2 denotes the dimensionless gamma ray photon energy, cr(x) ~ jqCf^x^/x 
for x > Xth = 2/[x (l — cosd)] is the pair creation cross section above threshold energy 
^thj °t the Thomson cross section and iV 7 (x ) = L a / (47irl xt m e c 3 x ) the external target 
photon density of the source with monochromatic luminosity L Q at target energy x Q . For 
L Q = 10 46 L 46 erg/s and x > x t h this yields 

r 77 (x) = 4 • 10 15 L 46 (rext/cm)- 1 x th x- 1 x- 1 (3.2) 

so that at threshold r 77 = 4 ■ 10 15 L^(r ext / cm)^ 1 x^ 1 . Assuming that the infrared emission 
comes from dust at r ext = 300 pc we obtain that r 77 > 1 for x > 5 • 10 5 corresponding to 
E 1 > 100 GeV. Since the produced pairs would isotropize rapidly, the reemitted radiation 
in the m.f. of the host galaxy does not remember the beam, so that the apparent luminosity 
is reduced to the comoving frame luminosity which is down by a large kinematic factor. 

One can further ask, wether the flattening of the predicted cascade spectrum above 
10 GeV (cf. Fig. [TJ) could be destroyed by photon-photon absorption as well as the 10 TeV 
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bump. The absorbing thermal optical/X-ray photons presumably come from behind the 
jet (from the apex) and therefore the angle 9 should have the value obtained for the angle 
between jet axis and observer, viz. 8 = 7°. With x denoting the observed gamma ray energy, 
the resonant target photon energy is x — 1.4 ■ (7 keV). The optical depth at threshold 
for L 46 (7keV) = 0.1 is given by r 77 (r c + r b = 8pc, 6 = 7°) ~ 6 • 1CT 4 ; clearly insufficient to 
steepen the spectrum. 

IV. CONCLUSIONS 

Radio jets emerge out of the vicinity of an accreting compact object, presumably a 
black hole, where a rotating magnetosphere tied to the accreting gas feeds a considerable 
fraction of the infalling matter into the outflow, thereby transporting angular momentum 
outwards. The magnetic fields allow for collimated flow solutions passing through bulk 
Lorentz factors of 5-20 beyond the lightcylinder at roughly 10rQ — 1.5 • 10 15 mg cm, where 
m 9 denotes the mass of the black hole in units of 10 9 m o . If the initial acceleration of 
plasma would be much more efficient, i.e. 71 ~ 10 3 , the ambient thermal UV and soft X- 
ray photons from the accreted gas excert a 'Compton drag' on the particles in the outflow 
slowing it down to terminal Lorentz factors m 10, thus dissipating most of its energy into 
Compton-scattered radiation [[K|]. However, this produces maximum photon energies well 
below 100 MeV - in apparent contradiction with the observed fluxes above GeV from many 
extragalactic radio sources. Moreover, if radio jets consist of ordinary hydrogen plasma, 
radio-loud objects like 3C 273 have a very powerful jet even many kiloparsecs away from 
the AGN and therefore it seems unlikely that they have dissipated much of their power in 
the central parsec already. 

After some inital collimation by toroidal fields the jet expands freely and super- 
Alfvenically into the surrounding medium. However, free expansion tends to bring the 
internal pressure of the jet within a few scale heights above the Alfven point into equilibrium 
with the external pressure. The flow tries to realise an equilibrium state by reconfinement 
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shocks ||, which let the internal pressure vary in a zig-zag fashion about the external pres- 
sure. Thus, the external gas determines the structure of the jet and the amount of radiative 
dissipation, which establishes itself through shock acceleration at the reconfmement edges. 
Another important source of shocks could be mass ejection at the base of the jet - much 
like in the solar wind. Thus, when the ejection sites are tied to a rotating accretion disk, 
a helical shock pattern results, which seems to be required by VLBI observations of many 
radio jets and HST observations of the M87 jet. 

Synchrotron emission of accelerated electrons emerges mostly in the radio to optical 
regime, with the submm/infrared/optical emission coming from the innermost nuclear jet 
(r < Tb ~ 4 pc) and the radio emission coming from further out. Note, however, that the 
propagation of the jet through the central parsec (out to the BLR) is assumed to be essen- 
tially non-dissipative, i.e. r c = 0[lpc]. During this collimated stage the initial Poynting 
flux is converted into bulk kinetic flux, which is then the reservoir for particle acceleration at 
shocks farther out. Higher frequencies than optical are difficult to achieve from synchrotron 
emission by the primary accelerated electrons, because the maximum energies are strongly 
constrained by energy losses ||. However, the protons accelerated in the jet reach energies 
so high (7 Pimax ~ 10 10 ) that saturation effects due to the finite size of the shocks become 
important, simply because their cooling time scale, which is dominated by secondary parti- 
cle production for a wide parameter range, is very long compared to the cooling time scale 
of electrons. Thermal photons from outside of the jet can additionally damp the protons 
in the jet. The power induced by the cooling ultra-high-energy protons at comoving frame 
energies up to 10 10 GeV is rapidly reprocessed by an electromagnetic shower and emerges 
in the X- to gamma ray regime. The infrared to optical emission from the throat of the jet 
serves as the scattering agent for the cascade determining the photon energy e 7 w 10 TeV 
where the jet becomes optically thick with respect to pair creation in the observer's frame. 

When the initial jet collimation terminates very far away from the central engine or 
when the thermal emission from the central engine is inherently weak, then the emerging 
cascade spectrum is much like as it is observed by GINGA, COMPTEL and EGRET for 
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3C 279: ax ~ 0.6 with steepening at a few MeV, so that a 7 1.0. On the other hand, a 
quasar like 3C 273 harbouring a "miniblazar" , but otherwise emitting predominantly thermal 
UV and IR emission, is very likely to show signs of interactions of these photons with the 
accelerated particles inside the nuclear jet. The interactions reduce the maximum energies 
of the particles they can obtain by diffusive acceleration limited by energy losses. As shown, 
this reduction moves the injection energy of secondaries at the top of the cascade closer to 
the energy where the jet becomes optically thick with respect to pair creation. Therefore 
3C 273 developes a different cascade spectrum with two bumps at a few MeV and TeV, 
respectively. The latter is absorbed by infrared photons from dust clouds within the host 
galaxy, whereas the MeV bump has power law wings on both sides in energy space matching 
the observed X-ray and gamma ray power-laws consistent with GINGA, COMPTEL and 
EGRET observations. Of course, this does not rule out the existence of additional emission 
components in the X-ray band. 

For 3C 273 flattening of the spectrum at a few GeV is expected, but observations do not 
seem to indicate this - a detail deserving further measurements. It is at the present stage 
of investigation unclear, what the variability properties should be like. Variability must 
be very complex, because we are actually dealing with a source extending over scales from 
fractions of a parsec to parsecs (transverse size of the jet or downstream emission regions 
of shocks). However, the variability time scale of the submm-nearinfrared target photons 
t\r — r b,_i_(l + z)/DjC ~ 180 hours) should at least be contained in the X- and gamma ray 
autocorrelation functions, the amplitude should be significantly larger because of cascade 
reprocessing. An outburst in the power of the jet itself should manifest itself first by a 
flaring sub mm/infrared/optical continuum (the target) followed very shortly after by the 
gamma rays and then much later by a corresponding rise of the radio flux. Since the radio 
flux peaks at its self-absorption frequency, the delay at a given radio frequency determines 
the scale length of the radio emitting part of the jet, viz. At ~ (r — r b )(l + z)j ' D s c with 
r = 4(z/ obs /z/ b ) _1 pc and v h — 2 • 10 10 Hz. 

As a corollary of the interpretation of the gamma rays from 3C 273 being of hadronic 
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origin it follows that 3C 273 should also be a source, which provides roughly 1% of a) the 
flux of cosmic rays at 10 19 eV at Earth and b) the atmospheric background neutrino flux at 
10 TeV in a solid angle of 1° x 1°. To come to these conclusions consider the following: 

(i) In contrast to the electromagnetic flux the neutron as well as the neutrino flux density 
spectrum from photomeson production are flat (the luminosity increases proportional with 
the energy), because the UHE protons in the jet scatter on a polychromatic target, i.e. an 
inverse power law, so that the number of target photons increases rapidly with increasing 
energy. The electromagnetic cascade washes out the injection energy of the photons and 
redistributes the power below the energy above which the jet is optically thick with respect 
to pair creation. 

(ii) Hadronic proton-photon interactions produce neutrons by isospin flip. The apparent 
neutron luminosity L n « L n + is given by L n « -^L 1 (Ref. |1| where pion and pair production 
are taken into account) and the maximum observed energy is E n = ^D- } E PjTDax / (1+z). Hence 
one obtains using fj F v dv = 3 ■ 10 13 Jy Hz 

F n _ p = ^£ n -Lx [ F^v ~ 7 • 1CT 19 cm-V 1 (4.1) 

while the observed flux at 10 19 eV is roughly 6-10~ 17 cm -2 s _1 [|l^]. The neutrons are not con- 
fined to the magnetic field in the jet and because of the /3-decay length lp = 100[7 n /10 10 ] kpc 
they can escape the host galaxy without adiabatic losses (cf. original papers in Ref. flOp. 
Energy losses on the cosmic microwave background further reduce the flux above 10 19 eV. 
Moreover, diffusion of the protons due to random magnetic fields in the surroundings of 
3C 273 could smear the cosmic ray beam, thus reducing the cosmic ray flux by a factor 
Dj~ 3 ~ 3 • 10~ 3 and the particle energy by D^ 1 = 1/7. Such smearing would, however, 
increase the total number of radio sources contributing to the entire extragalactic flux, since 
then also non-beamed (not core dominated, steep spectrum) sources would contribute (cf. 
Ref. |T7| for a detailed model with direct proton escape from isotropic hot spots at the tip 
of powerful jets). . 

(iii) The neutrino luminosity for all species is given by L u « fs^-y an d E Utmax = 
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2^-Dji?p >max /(l + z), see Ref. [|T^]. This yields the flux 



F„(p7) = ^E-^ jT F„di/ ~ 5 ■ 10- 18 cm^s- 1 (4.2) 



important at Fly's Eye energies |Tj|. The neutrino flux from pp-interactions (which domi- 



nates at energies in the 10 TeV range, see Ref. ||) is given by 
^ W - TlTi^ Fu(pr)(lA-) ^2-10- 5 m 9 F,(p 7 )f inin i l„,^ 1 (4.3) 



rvlnF IE ■ v U \E I ~ v 7 V 10 10 GeV 

reaching 10~ 16 mg cm -2 s _1 at 10 TeV. Thus, the neutrino flux from the jet of 3C 273 
is far below the background neutrino flux from AGN as predicted by several au- 



thors (see Refs. |19| , p0| ). However, the effective intensity for a solid angle of 1° x 1° 
reaches 3 • 10~ 17 m 9 cm" 2 s _1 sr _1 GeV -1 at 10 TeV for 3C 273 alone, compared to 
~ 10" 14 cm -2 s _1 sr _1 GeV -1 for the hypothetical AGN background (which itself becomes 
comparable with the atmospheric background at 10 TeV). Considering that the number 
of similar flat-spectrum radio sources is very large [El| detection of the entire population 



of these powerful cosmic accelerators seems possible by refined techniques using horizontal 



atmospheric showers [23] and by taking advantage from the fact that the candidate source 



positions are known to any desired accuracy. 
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FIGURES 

FIG. 1. The broadband spectrum of 3C 273 from submm wavelengths up to hard gamma rays 
decomposed into a "blazar" component from a nuclear jet and thermal emission from accreted 
gas. Data for the submm-optical blazar component are taken from Ref. ||, who suggested the 
decomposition based upon polarization and variability, for the X-ray range from Ref. |p3| , for the 
hard X-ray/soft gamma ray range from Ref. |24|] and, finally, for the gamma ray range from Ref. 
jlj]. Several attempts to detect 3C 273 with Cherenkov telescopes have failed p5j| . The synthetic 
spectrum shown comprises the emission from shock accelerated electrons (submm-optical) and 
protons (X-gamma rays). The emission mechanism is optically thin synchrotron emission from 
"primary" (=shock accelerated) electrons and from pairs produced by an unsaturated synchrotron 
cascade initiated by the protons cooling in the soft synchrotron photon soup. The TeV emission 
is absorbed by infrared photons from dust within the host galaxy. The model predicts a flattening 
of the spectrum between 10 GeV and a few hundred GeV. 
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TABLES 

TABLE I. Physical conditions in the compact jet of 3C 273. Basic parameters for the cal- 
culation of physical conditions are: redshift z, Lorentz factor of jet 7j, angle to the l.o.s. 8, 
proton/electron energy ratio r/ = Lj/L^ r / a , equipartition constant k c ~ u re [/uB, luminosity of jet 
in m.f. and rel. particles L44 in units of 10 44 erg/s, observed opening angle of jet $ Q b and v c 
comoving frame cutoff frequency of the synchrotron spectrum from the primary electrons. Inferred 
quantities (consistent with observations) are the observed apparent speed /3 Q b, the Doppler factor 
of the jet Dj = [7j(l - $ cos6>)]~\ the primary electron break frequency 14,, the photon/magnetic 
energy ratio a s due to the radiation produced in the jet and due to an external luminosity of 
10 47 erg/s, resp., the projected length of the radiative jet rb,ob where most of the emission with 
v > 1% comes from, the deprojected length r^, the (unresolved) transverse size of the jet rj_, the 
target photon compactness in the comoving frame I, the magnetic field at r^, the electron and 
proton break Lorentz factors and 7 p b and finally the observed photon energy E* where the 
jet becomes optically thick with respect to pair creation. 
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